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Abstract 

The masses of excited heavy mesons are studied with sum rules in the 
heavy quark effective theory. A set of interpolating currents creating (anni- 
hilating) excited heavy mesons with arbitrary spin and parity are proposed 
and their properties are discussed. Numerical results at the leading order of 
the 0(l/mg) expansion are obtained for the lowest doublets (0 + ,l + ) and 
(l+,2+). 
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Properties of excited heavy mesons composed of a heavy quark and a light anti-quark 
are of interest for several reasons. Some of these excited states have been observed in 
experiments. They will be objects of further study in future experiments with B-factories. 
In particular, they are useful for tagging in CP experiments 0. Theoretically, the rela- 
tive simplicity of the dynamical condition makes it of interest for exploring the internal 
dynamics of systems containing a light quark. 

The I/tuq expansions are useful for such systems (mq is the heavy quark mass). The 
heavy quark effective theory (HQET) |2|] is a powerful tool. However, for obtaining some 
detailed predictions one needs to combine it with some non-perturbative methods. The 
spectra and decay widths of heavy meson excited states have been studied with the 1/tuq 
expansion in the relativistic Bethe-Salpeter equations in [3,4]. They can also be studied 
with the QCD sum rules in HQET which have been used for ground states of heavy 
mesons [5, 6]. 

As a first step, in the present article we shall use the QCD sum rules in HQET to 
obtain the masses of lowest heavy meson doublets (0 + , 1 + ) and (l + ,2 + ). These states 
have been studied with the sum rules in the full QCD in the literature @. We shall see 
in the following that using the HQET has some advantage for such problems. 

In || it was shown that if the interaction kernel contains only 7 M and I on the heavy 
quark line, solutions of the Bethe-Salpeter equation for heavy mesons of arbitrary spin 
and definite parity consist of two series of degenerate spin doublets in the infinity quark 
mass limit tjiq — > oo. The wave functions XjPii where j, P are spin and parity and i — 1, 2 
for these two series of states, have the following most general form 
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In the above formulas, is the meson mass, x is the coordinate of the heavy quark, u 
is the coordinate of the light anti-quark relative to the heavy quark, v is the velocity of 
the heavy meson, u t = v ■ u, u til = u M - up^, 7 iM = ln~ Av> Vai-a, is the transverse, 
symmetric and traceless polarization tensor. Finally, and \&i are invariant functions of 
ul and U{. 

The assumptions used to obtain these results can be slightly generalized compared to 
that explicitly written down in M. We only need to assume that in the uiq — > 00 limit 
the kernel is of the form 

G(P, p, p') = I ® r U S (P, p, p') + 7 M ® r M U v (P, p, p') (3) 



where P is the total momentum, p and p' are the relative momenta of the quarks in the 
initial and final state respectively, V and r M are arbitrary functions of momenta and 7 
matrices for the light quark which transform as a scalar and a vector respectively, Us 
and LV are arbitrary invariant functions of momenta. Instead of the free propagator, 
the general form can be used for quark propagators appearing in the B-S equation. The 
calculations leading to above results only need to be slightly changed for this general case. 
The restriction to / and 7^ in (|3]) for the heavy quark is equivalent to that the interaction 
is independent of the heavy quark spin in the limit vtlq — > 00. This is satisfied by QCD.Q 
Besides this, the interaction kernel is most general. The B-S wave function is not gauge 
invariant so that functions $j and are gauge dependent. However, the general form of 
(1) and (0) is the same for any gauge. Physically, each doublet has a definite value for the 
total spin of the light component. Except for the states 0~ and + , for each state in the 
series (1) there is a state with the same j, P in the series (0). They correspond to different 

values of ji = j + - and j — -. In the following we shall switch between the indices i 
and ji. This pattern is in agreement with the heavy quark symmetry. The above results 
are model independent consequences of the heavy quark symmetry. The B-S equation, 



1 In QCD the elementary interaction vertex of the quark with the gluon is 7^, the heavy quark 

1 1+ i> ■ 1 ,. . tt . 1+ i> 1+ i> 1+ f , . 

propagator is m the uiq — ► 00 limit. Using 7^ = , the contribution 

i) ' k 2 2 2 2 

of any Feynman diagram to the B-S kernel can be reduced to the form (^) after integrating out 

the internal momentum of the diagram. 
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being Lorentz covariant, only serves as a convenient tool for finding the form of the wave 
function of the state with definite values for the angular momenta j, ji and the parity 
P. For related results, see PJlOj. The authors of these latter works used directly the 



transformation properties of the wave function of a composite system with the angular 



momenta j and ji . \W\ contains results similar to ours, though the author did not deal 
exactly with the B-S wave function. This again implies that our results on the forms of 
the B-S wave functions in the tjiq — > oo limit are model independent. 

The relations between B-S wave functions and the matrix elements of the currents 
make the above results useful for constructing the appropriate interpolating currents and 
deriving their properties. The wave function XjPi i n (1) an d (0) can be written in the 
form 

Xj, P ,i = ^•^r 7a ,.. a T^(A) Uiulm) +il-f 2 faui) \ . (4) 



We propose to use 



(■<•) = hwrfep^vMx) , (5a) 



or 



•/•;:,■" (.'■! = hAr)\j r ;r ] (n,n ») v^) , (5b) 

as the interpolating current which creates the corresponding heavy meson, where 
h v ( x ) = e imBV - x ^-^-Q(x) is the heavy quark field in HQET and Q(x) is that in full QCD, 

d 

T> Xt is the covariant derivative — igAAx) and 

dx t 

V UK n ] (V,..) = Symmetrize {i—^J - ^g^g^T^"^ (P Xt )} 
with g? ia2 = g aia2 - v ai v a \ 

The currents hermitian conjugate to (Ba|) and ( 55[ ) 

•/;:;:7' i-'-) = q{x)Y { ^p\v Xt )h v(x) , (6a) 
= V x ^tpp\v xt )h v{x) , (6b) 

where f = 7°r^7°, annihilate the same meson. The currents defined in (5) and (6) have 
nice properties. Let \j,P,i) be the heavy meson state with the quantum numbers j, P, i 
in the niq — > oo limit. We have 
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where fjPi = fpj t has the same value for the two states in the same doublet. The formula 
(0) can be verified as the following. Without loss of generality we can choose a gauge 
satisfying A^(x ) = 0. In such a gauge, it follows from the definition of the B-S wave 
function Xj,P,i(. x > u ) = (Q(x)q(x + u)) \j, P, i) that for a current in (|6a|) 
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where A = — tuq. Using the forms of V in (1) and (2) and the relations of the type 
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where S denotes symmetrization and subtracting the trace terms in the sets (a% ■ ■ -aj) 
and (Pi ■ ■ ■ Pf), (0) can be obtained from (§). Only the term fi contributes to fpj r If 
instead of J we use J' in fl6b|), the calculation is almost the same, except that now only 
the fi term contributes to fpj r 

Using the equation of motion for the light quark ij)q = and that for the heavy quark 
in HQET iv ■ Dh v = 0, one can obtain the relation 



(10) 



between the coupling constants for the current with the extra factor ij) and that without 
it. 

Similarly, calculated with the quark and gluon fields in the mg — > oo limit the corre- 
lation functions of Jj,p,i and Jl P , v satisfy 

i <o|t (./; , ;;v'' (,)./;/;;;;;'((.)) | > = %,<wM-i) j s ■ ■ 

x J dt5(x-vt)U Pjjl (x) , (11) 
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where Hpj t is the same function for the two states in the same doublet. For verifying (11) 
we use the Schwinger gauge 



= 0. 
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The heavy quark propagator in HQET is 



i (0\T (h v (x)h v (0)) |0) 



1+ jt 



dt5 i (x - vt) . 



(13) 



The leading order Lagrangian in HQET is Cq = ih v v ■ T>h v . Therefore in gauge (|T2"D there 
is no gluon line terminating on the heavy quark line. We thus obtain 
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Using the relation 
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in the gauge (0), one can rewrite (|T4|) in the form 

i (0\T (jj,p,i( x ) J j',P'A )) 1°) = ~ l j dtS4 ( x ~~ w*)*r{r iW / 
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From Lorentz and translation invariance for gauge invariant quantities, S(y, x) has the 
general form 



S(y,x)^-^-=\Sx((y-x)^,yi-xi)+i ^^((y-x^^i-xiy ' 



(17) 



Substituting ( IPTf ) into dl5|) brings it into a form almost identical to the r.h.s. of @. The 
formula (|ll]) then follows after similar calculations. 
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(^D and ([TTf) imply that two currents in the sets (|5a|) and ( |5bl) with not identical values 
of j, P, i never mix in the mg — > oo limit either in the heavy meson or in the quark-gluon 
level and heavy quark symmetry is explicit with these currents. This verifies that they 
are the appropriate interpolating currents for heavy meson states with definite j, P and 
the light quark momentum ji which are conserved in QCD in itlq — > oo limit. 

The properties for currents mentioned above are important for applications to QCD 
sum rules for excited heavy mesons. If one use other currents or study the sum rules 
in full QCD, there are in general contributions from two nearby poles corresponding to 
states of the same j, P. Their contributions may not be separated correctly. Furthermore, 
the mixing of such two states can be calculated within our formalism by introducing the 
OiXjvTLQ) terms in the Lagrangian in HQET. 

For obtaining QCD sum rules in the HQET we study the correlator 

n^r^^W = ijd 4 xe ik -*{o\T (7;;;;7 1 '(.r)./; ;;: '(())) |o> . (is) 

From (|^), it has the following form at the leading order 

K p T A '"^ (*) = & 9?* ■ ■ ■ 9? Pi IWw) , (19) 

where u = 2v ■ k and Uj p ti (u) = Upj^u) is independent of j. 

We shall confine us to the leading order terms in this article. The contribution to 
Upj^uj) of the pole term on the hadron side is 

On the quark-gluon side of the sum rule, Upj^uj) is calculated with the leading order 
Lagrangian Co in the HQET. The Feynman diagrams contributing to ([18]) are shown in 
Fig 1. For the sake of simplicity we consider the doublets (0 + , 1 + ) and (l + ,2 + ) in the 
following. 

According to (|5a| ) and (f[5|), there are two possible choices for currents for the doublet 
(0+1+), either 

4+,a = ^ M > ( 2 °) 
J\% 2 = -L h vl 5 ^ q , (21) 

or 
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Jo]+,2 = ^=K{-i) Ptq, (22) 
Ji]+,2 = ^h v l 5 l?H) Ptq- (23) 

Similarly, there are two possible choices for the currents for the doublet (1 + , 2 + ). One is 

= \/f WH) - k Q A) g , (24) 



. 4 v ' V 1 3 
4,+T = 1/f {itV? + 1?V? - 2 9 r 2 Pt) q ■ (25) 



Another choice is obtained by adding a factor —i J) t to fl24|) and (p5|). Note that, without 
the last term in the bracket in (23) the current would couple also to the 1 + state in the 



doublet (0 + , 1 + ) even in the limit of infinite mq. 

Usually, the currents with the least number of derivatives are used in the QCD sum 
rule approach. The sum rules with them have better convergence in the high energy region 
and often have better stability. However, there is a motivation for using the currents (f2~2|), 
(P5|) for the doublet (0 + , 1 + ). In the non-relativistic quark model, which usually gives 
correct ordering of energy levels of hadron states, the doublets (0 + , 1 + ) and (l + ,2 + ) are 
orbital p-wave states which correspond to one derivative in the space wave functions. In 
fact, as mentioned above, the coupling constant for the currents is proportional 



to ^2 in (2) and that for the currents (24), (25) is proportional to $i in (1) which are 
the large components in the non-relativistic approximation. Therefore, we shall consider 
both the currents (|20j), (gTJ) and (|22|), fl2|) for the doublet (0 + , 1+). 



First we consider the doublet (0 + , 1 + ). For the correlator 

n , +i2 = i J d'xe^(0\T (j' , + ,2(z)4,+,2(0)) |0> , (26) 

where J'0,+,2 is the current conjugate to fl22|), the sum rules after the Borel transformation 
is found to be 

f 2 e~ 2A/T = ^ jf wV^cfc, - ^ml(qq) , (27) 



• 5 — ' T du 



3 f u)c 

^3T2 / U " e ~" dUJ 

In the above sum rule we have confined us to terms of lowest order in perturbation and 
operators of dimension less than six. 
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The correlator of the 1 + current J'j.,+,2 (H) and its conjugate has the form 

rr^ i2 = < / d 4 xe^(o|T (j? >+>a (x) J't + M) 1°) = n w ( w ) . (29) 

We find for n 1+i2 (u;) the identical sum rule as that for n 0i+i 2(<^), as it should be. 

The corresponding formula when the current Jo,+,2 and Ji,+,2 in (|20|) and (^) are used 
instead of «/'o,+,2 and J'1,+,2 are the following 

/ 2 e" 2A/T = ^ jT ^e-^du + i(gg) - ^(g,) , (30) 

A = 3 ^ ; ^ (31) 

— J o u 2 e u du + (qq) - -^m 2 (qq) 

The correlators for the (1 + , 2 + ) doublet have the form 

K^k) = -gf n + , jl=3/2 (u) , (32) 

nSffi AA (*) = ^ (<? t Ql/ V 2/32 + flT*^ - ^r Q2 ^ lA ) n +>ii = 3 /2H . (33) 



When the currents (|24|) and ( pq) are used the sum rules for this doublet are found to be 

Here the ( — G 2 ) term comes from the diagram (c) in Fig. 1 which vanishes for the doublet 
(0+1+). 



A= ; 25 r J0 ~~2 ~ ~-, ~ • 05) 



From (P8|) , (pl|) and (p5[) we obtain curves shown in Fig. 2a,b,c. For the first two cases 
stability windows exist when u c lies in the interval 2.4 to 2.9 GeV; for the last one uj c lies 
in the interval 2.7 to 3.2 GeV. Imposing usual criterium for the upper and lower bounds 
of the Borel transform parameter T we found approximately the same stability window 
0.7-0.9 GeV. The results for A for three cases are the following 

A(l/2 + ) = 1.05 ±0.05 GeV, (36) 



for the doublet (0 + , 1 + ) when the currents (|20[ ), fl2lD without the derivative is used 
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A(l/2 + ) = 0.90 ±0.10 GeV, (37) 
for the same doublet when the currents (23), (|2"B|) with the derivative is used. 



A(3/2 + ) = 0.95 ±0.10 GeV, (38) 

for the doublet (l + ,2+). 

The predicted mass value for the doublet (0 + , 1 + ) in ( |36|) is somewhat different from 
that in (|37|). The value in ( j36|) is in agreement with the result of [7] in which the same 
currents as (|^) and (pTj) are used. (|36|) predicts a mass value of the doublet (0 + , 1 + ) 



slightly higher than the mass value of the doublet (1 + , 2 + ) predicted in (^). On the other 
hand, ([37]) gives a mass value for the doublet (0 + , 1 + ) approximately equal or slightly 
lower than that for the doublet (1 + , 2 + ) in (p8|). Though (36) and (37) can be considered 
as consistent within the uncertainties, (37) indicates that the QCD sum rule approach can 
incorporate a mass value for (0 + , 1 + ) lower than that for (l + ,2 + ). This is in agreement 
with quark model predictions in [11] and [3]. By contrast, we have also derived the sum 
rule for the ground state doublet (0~, 1~) with the currents with the extra derivative —i'p. 
We found that there is no acceptable stability window for this sum rule. 

Finally, we would like to address the following points of this paper. The results from 
our previous work on B-S wave functions are used only for finding the appropriate currents. 
Once these currents are given, the calculations in QCD sum rules are independent of the 
consideration of the B-S equation. That the currents used are appropriate was checked 
by proving (11) with direct calculations in terms of quark-gluon fields in HQET without 
referring to the B-S wave functions. Assuming duality and inserting intermediate heavy 
meson states between J and J* in (11), this equation also strongly suggests equation (7) 
which is the direct statement that the currents create pure states with definite values of 
j, P and j t . 
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Figure Captions 

Fig. 1. Feynman diagrams contributing to the sum rule. There is no 
Fig (c) when the currents (20) (21) are used. 

Fig. 2. Dependence of A on the Borel parameter T for different values 
of the continuum threshold uj c . (a)(b) for (0 + , 1 + ) doublet with the 
interpolating currents (|22|) ( p3| ) and (|2^) ( |2T| ) used respectively. From 
top to bottom the curves correspond to (a) uj c = 2.8, 2.6, 2.4 GeV and 
(b) uj c = 2.9, 2.65, 2.4 GeV. (C) for doublet (1+, 2+). From top to bottom 
the curves correspond to uj c = 3.1,2.9,2.7 GeV respectively. 



11 



Fig 1 
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